Interfacial instability and DNA fork reversal by repair proteins 
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A repair protein like RecG moves the stalled replication fork in the direction from the zipped to 
the unzipped state of DNA. It is proposed here that a softening of the zipped-unzipped interface at 
the fork results in the front propagating towards the unzipped side. In this scenario, an ordinary 
helicase destabilizes the zipped state locally near the interface and the fork propagates towards the 
zipped side. The softening of the interface can be produced by the aromatic interaction, predicted 
from crystal structure, between RecG and the nascent broken base pairs at the Y-fork. A numerical 
analysis of the model also reveals the possibility of a stop and go type motion. 



I. INTRODUCTION 



The semi-conservative replication of DNA requires its 
unzipping by helicases, and synthesis of new strands over 
the opened parent strands by dna-polymerases (DNA- 
polymerase III) with high processivityP, The basic 
features of rephcation arc as follows. The two strands of 
DNA run in opposite directions. The polymerase with 
unidirectional 3'-5' motion can then follow the helicase 
on one strand (the leading strand) only, but not on the 
other (the lagging strand). Though the two polymerases 
are identical, the lagging strand polymerase performs a 
more complicated job, e.g. (i) it makes the new strand 
in pieces (Okazaki fragments), (ii) it has to shift repeat- 
edly to newer positions closer to the moving helicase, (iii) 
it restarts the polymerization process anew after every 
shift, and, then (iv) the small Okazaki fragments are to 
be joined by dna ligase. In addition primase is required 
for initiation of the Okazaki fragments, sliding clamp pro- 
teins to tether the polymerases to DNA, to name a few 
more. Moreover, the fidelity of replication requires addi- 
tional repair or proof-reading capabilities which need to 
become functional as and when needed. Processivity in 
this context is defined as the number of base pairs added 
at a stretch each time the replication machinery binds to 
DNA (for E. Coh, 4.6 x 10^ base pair long DNA is ulti- 
mately replicated in 40mins) . It is understood that most 
of the polymeric molecules of the replication machinery 
can work independently in vitro but all of these need to 
work in tandem over a long time and long distance (along 
the DNA) for the processivity observed during replica- 
tion. 

The controversy with the replication dynamics can be 
summarized in the following way [11, In one class of 
models, the replication machineries move on the DNA. 
See Fig. [1^. One needs tight coupling between the lag- 
ging strand polymerase and the helicase, though they 
are independent, to facilitate repeated recognition of the 
newly unzipped region of DNA. Not only the polymerase 
and helicase, others like the primase, ligase also need to 



get correlated on the lagging strand. The antithesis is 
the proposal of a replisome - a complex of all the objects 
staying together with the DNA looping through the com- 
plex in a very particular nianner|i3|. See Fig. [T}d. Since 
the DNA goes through the replisome, there is expected to 
be a depletion layer of nutrients surrounding a replisome 
and it needs to be supplemented by a current towards it. 
More perplexing is that, if a complex can exist for a long 
enough time during replication, then why it is so elusive 
in vitro and what is responsible for the "bound state" in 
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FIG. 1: Schematic diagram of the replication process. Only 
the helicase and the polymerases are shown with arrows in- 
dicating directions of motion. The solid lines are the parent 
DNA strands with 3', 5' ends marked. The dash-dot lines 
represent the newly synthesized DNA. In (a), the members of 
the rephcation machinery move on the DNA individually. The 
small fragments on the lagging strand are the Okazaki frag- 
ments. The replisome model is shown in (b). The small thick 
lined box is the clamp loader that loads the lagging strand 
DNA on the replisome. In both models, the whole machinery 
works near the Y-fork. 



Single molecular experiments have revealed correla- 
tions between synthesis by T7 DNA polymerase and T7 
helicase activity on a dsDNA, e.g. efficient duplex DNA 
synthesis require combined action of polymerase and he- 
licase but no direct or specific interaction between the 
polymerase and helicase play any roleQ. For bacte- 
riophage T4, the assembly of the polymerase and the 
clamp-loader with DNA has been found to be highly dy- 
namic and well coordinated Q. The "replisome" may it- 
self be dynamic in nature with not just two but even 
three polymerases, and others like repair factors can be 
dynamically attached to itQ- In case the replication 
process were tightly controlled, hard-coded in the func- 
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tioning of the molecules, then that would be reflected in 
the statistics of the lengths of the Okazaki fragments on 
the lagging strand as well. The probability distribution 
of the lengths of these fragments for T4 and T7 have 
been found to be very broad and highly non-gaussian, 
compared to very sharp distributions expected from in- 
dividual primase controlled models [1]. Purely stochas- 
tic models have been proposed to study correlations in 
Okazaki fragment distributions @. These newer experi- 
mental methods and analysis, apart from elucidating the 
nature of the complexity of replication, are also bringing 
out the importance of dynamics and correlations (also 
called coordination) of the protein complexes and DNA 
strands in particular near the fork. 
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FIG. 2: Schematic diagram of the repair process by RecG 
(shaded ellipse and two circles) at the Y-fork. (a) RecG at 
the stalled Y-fork with a longer lagging strand pair, (b) Fork 
reversal and chicken-foot configuration (Holliday junction). 
The arrow denotes the direction of motion of the fork. The 
four arm junction is created by the fork reversal and the con- 
tinuation of the leading strand pair by making use of the pair 
on the lagging strand. 

An important aspect of the replication process is the 
repair mechanism also involving helicases. A helicase is 
a motor protein that generally facilitates in opening the 
DNA and leads the replication machineries. Such a he- 
licase moves from the unzipped (i.e., the open side) to 
the closed side. There are actually many helicases, e.g. 
almost 20 in E. Coli (though not clear why so many), 
but they generically are known to have a part (called 
the wedge domain) that maintains the two strands of the 
DNA at a distance much larger than the base pair dis- 
tance. The motor action (e.g. for dnaBQ) and often 
additional pulling action (e.g. PcrAp^) carry the open- 
ing process on. It stops in case there is a nick or a lesion 
on the leading strand. The replication process stalls, the 
whole replication assembly disperses, and a repair pro- 
cess takes over. Often, RecG is the helicase that starts 
the repair job by performing a fork reversal Based 
on the crystal structure, it was suggested that the wedge 
domain maintains the separation but in addition there 
arc two domains that sit on the zipped side and pushes 
the DNA towards the unzipped side thereby scooping 
out the DNA on the lagging strand[ll|- See Fig. [H This 
continues until the nascent end on the leading strand is 
reached, forming a chickenfoot configuration. The job of 
RecG is over, the repair process is completed by other 
objects via a Holliday junction and then the replication 
restarts. In this model, RecG has two jobs, zipping the 
parent strands and unzipping the nascent duplex on the 



lagging strand. The activities of a pre-asscmbled RecG- 
DNA on (i) a three way DNA (consisting of the Y-fork 
with a nascent duplex on the lagging strand) and (ii) a 
two-way DNA (consisting of a pure Y-fork) are found to 
be the same, if initial transients are ignored p^. 

We therefore take this view that the main job in the 
repair process of RecG is the zipping of the DNA or the 
fork reversal. 

Our purpose in this paper is to show how a coupling of 
the helicase and the replication fork (Y-fork) can be used 
to formulate the backward mobility of RccG-like repair 
proteins, in addition to the forward motion of general 
helicases. Our hypothesis is that the Y-fork is an essen- 
tial element in the whole replication pro cess, and it is at 
the heart of the correlated dynamics [IJ]. This hypothe- 
sis is based on the observation that in both the biological 
models and experiments mentioned earlier, the common 
feature is the proximity of a Y-fork - the junction of an 
unzipped and a zipped DNA, and dynamics plays an im- 
portant role, may even be responsible for the formation 
of the replisomc. 

The first step to motivate this connection is to ana- 
lyze the helicase activity. Traditionally, DNA opening is 
considered a melting phenomenon where thermal noise, 
fluctuations in base pair breaking-joining, and polymer 
configurations play important roles A different mech- 
anism is the force induced unzipping transition at tem- 
peratures below the melting point[15i. ,16. ,17,. .18.]. In the 
unzipping scenario, a coexistence of the two phases sep- 
arated by a domain wall, at this first-order transition, is 
the Y-fork. The wedge domain of a helicase provides the 
constraint to maintain the DNA in a fixed distance en- 
semble, and, therefore, a Y-fork develops. The helicase 
activity is then tantamount to setting the Y-fork inter- 
face in motion. In this language, the replication problem 
can be recast as a velocity selection problem. The Y-fork 
dynamics is described by the internal dynamics of DNA, 
the helicase motion is controlled by the energy supply and 
its own dynamical mechanism, and similarly for other 
individuals. Thus, each has its own characteristic veloc- 
ity acting alone, and that's the velocity one measures in 
vitro [l^- Despite this diversity, these should have the 
same velocity when together during replication. Hence, 
the issue of velocity selection. Similarly, a repair process 
is not just a resealing of the Y-fork under DNA dynam- 
ics but has to be closely knit with the repair process. 
The specialty of RecG comes from the crystal structure 
which indicates a competition of the closed rings of the 
wedge domain with the nascent broken pairs of the DNA 
through aromatic interaction[ll|, That both the for- 
ward and the backward motions can be handled in the 
same framework lends credence to the basic hypothesis. 

The outline of the paper is as follows. Using a Landau 
type functional to describe the zipped-unzipped coexis- 
tence we formulate the propagating front equation when 
there is a local perturbation around the interface. The 
coupling of RecG with the interface is then introduced 
and the effective dynamics written down. This strong 
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coupling that can lead to the elastic or diffusive term 
to change sign, stabilized by a fourth derivative term, is 
then studied. This formulation is done in Sec II. A per- 
turbative approach is used to calculate, via a Goldstone 
mode, the velocity of propagation, the starting point be- 
ing the zero velocity coexisting phase situation. This is 
done in Sec III. A numerically exact solution of a dis- 
cretized version of the model is then presented in Sec IV. 
Apart from verification of the predictions of the pertur- 
bation theory, nonperturbative effects are also found. Sec 
V is the summary and conclusion. 



II. MODEL AND VELOCITY SELECTION 

A. Dynamics 

We start with the fact that a double-stranded DNA 
in equilibrium below its melting point can show phase 
coexistence with an interface separating the two phases, 
the zipped and the unzipped phases. A helicase (or its 
wedge domain mentioned in Sec I) provides the necessary 
constraint of a fixed distance ensemble to generate the 
coexistence. To open the DNA this interface needs to be 
translocated. 
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FIG. 3: (a) vs for /i = and h = 0.4. The perturba- 

tion by h makes the unbound state the stable phase, (b) a(fc) 
vs k for the D' < case. The minima at fc = ±fco suggest that 
the stable state is not the homogeneous state with = or 1, 
but a different state with a modulating wave- vector k = ko- 

For standard helicases like dnaB, the motor action just 
forces through the DNA so that as the location of the con- 
straint of the fixed distance ensemble shifts, the bound 
state region becomes metastable or unstable with respect 
to the unzipped state. For helicases like PcrA, there is an 
additional pull by a hand-like branch on the DNA strand 
near the Y-fork. This is also an example of an active 
mechanism to make the state near the Y-fork meta- or 
unstable but not in the bulk. With this in mind, let us 
formulate our model in the following way. The equilib- 
rium phase coexistence - Y-fork - by a static helicase, can 
be described by a Landau like free energy Fh{4') where 
(/) is a scaled distance. By choosing (/> = as the bound 
state and </> = 1 as the open state, with a slightly negative 
(j) representing an overtight state, we take 
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where D is an elastic constant and 
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The phase coexistence occurs at /i = 0. The coefficients 
are chosen, for simplicity, to have the extrema of in 
a symmetrical fashion. The perturbation h > makes 
the unzipped state the favoured phase. See Fig. [3^. The 
dynamics is described by 



with the boundary condition: 

0(-cx)) = 1, 0(+oo) = 0, 



(4) 



(5) 



and h{z, t) is non-zero only near the interface. In Eqs. [T] 
and m //i, represents the effective force. The bound- 
ary conditions ensure that we have open or unzipped 
strand on the left side while a closed terminal on the right 
side. The chain-length has been taken to be infinity. De- 
spite the resemblance, this is not the Fisher-Kolmogorov 
problem [23, [Ml because there is no bulk instability here. 

Where is the helicase in this set of equations? For a 
static case, ft, = 0, and there is no propagation. The 
boundary conditions impose a kink 4>{z) that goes from 
0=1 to = 0, but the location zq of the kink is 
arbitrary [22. This is a Goldstone-like mode; the inter- 
face or the kink separating the unzipped and the zipped 
phases can be placed anywhere along the chain because 
neither Eq. [5] cares about zqj nor there is any energy cost 
in shifting the interface. However, a static helicase fixes 
the position of the Y-fork and therefore it kills the Gold- 
stone mode by fixing zq where the Y-fork would be. The 
local instability is taken into account by h{z, t) which 
parametrizes the motor action and the active mechanism 
of the helicase. In short, the helicase at time t = is re- 
placed by two features: (i) the time independent bound- 
ary conditions and the location of the initial Y-fork rep- 
resent the wedge domain, and (ii) the passive or active 
process of the helicase is represented by the local insta- 
bility parameter h. 

The dynamic activity of a helicase makes zq a function 
of time. A propagating mode would imply 



(z,t) = $(z-rf), 



(6) 



■ so that the effective Landau-Hamiltonian can be written 



with c the velocity of the front. The helicase must move 
with the same velocity and the perturbation too, so that 

h{z,t) = H{{z - ct)/W), H{z) ^ for \z\ < 1. (7) 

This matching of velocity is the velocity selection men- 
tioned earlier. W is the width of the region over which the 
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helicase affects the fork. A simpler, more practical, choice 
would be an implicit definition h ^ Q for |(/) — < 5(j)yj. 

RecG is similar to other helicases so far as the exis- 
tence of a wedge domain is concerned, but its interaction 
is with the broken bonds at the interface (more like a 
surfactant). In other words, unlike the ordinary heli- 
cases that unzips, it does not produce any instability of 
any phase but rather interacts with the interface. With 
that in mind we introduce a new Gaussian variable y(z, t) 
representing the helicase and its interaction with the in- 
terface (for which d(j)/dz ^ 0) as 



^{0} 



dz 



Fo(0) + 2Ay- 



D 

'2\dzJ ' ' ~"" dz 

(8) 

with (K > 0), where the subscript indicates h = in 
Eq. [TJ i.e., there is no perturbation in the coexisting part 
of the free energy. 

Instead of coupled dynamics, we consider the effective 
dynamics of the Y-fork by integrating out the y field. 
Such an integration leads to a new effective Hamilto- 
nian of the type Eq. [21 with a reduced elastic constant 
D' = D — /K . For sufficiently large A, D' can become 
negative and for stability a higher order term is added. 
The effective Landau energy can then be taken as 
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The dynamics is given by 



dt 
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)(l-0), (10) 



where the helicase (perturbing) part is written explicitly 
as the hf( term with hf( ^ in a region near = 0. 
We use the fact that RecG operates near the zipped side. 
With negative D' , there is a preference for a modulated 
structure. To see the effect of D' < 0, we rewrite the 
gradient-dependent part of J-{4>} in Fourier modes as 



where 



(^) 



dk 

(27r)i/2^ 



dk a{k) (j)k(j>-k, (11) 



and a{k) = D'k^ + jk^ 



This favours, as shown in Fig. [SJd, a modulated state with 
wave- vector fcg = j/\D'\. Such a modulated structure 
would be a state with bubbles - but that does not hap- 
pen here because the perturbation is only local, meant 
to destabilize the interface. 

The difference between the two types of helicases is 
now seen in the respective /i-type terms. In Eq. [51 the 
/z-term affects the stability of the phase while in Eq. [TUl 
it affects the interface. 



III. PERTURBATIVE APPROACH 

We now treat the helicase effect as a small perturbation 
and show the change in the direction of velocity in the 
two cases in a first order perturbation theory. Let us take 
Eq. [2] for illustration. If there is no perturbation, then 
the free-energy, Eqs. HI suggest a static interface, 
i.e.; velocity c = 0. The profile with small perturbation 
can therefore be written as 



(j) = <l)o{z) + S<p{z — ct, t) 



(12) 



where the zeroth order solution is the static solution sat- 
isfying 

D^ + U<Po-]:){l-M^Q- (13) 
oz'^ 2 

For the boundary conditions of Eq. [H Eq. [T31 gives the 
static kink solution (j)n(z) that gives the z-depcndent pro- 
file near the interface |23|. The perturbed part satisfies to 
first order an inhomogeneous differential equation 



d 9^ 



S4' = Sf(4>o) = hMi-4'o) (14) 



where /'(0o) ~ df / d(j)\,f,=(i,g . This solution of this equa- 
tion can be written in terms of the Green function 12111 



dr 

tn J —OO 



d^Git,T,z,0 SfiM (15) 



Moreover, at first order, the velocity c is small and the 
perturbed part can be written as 



dz 



{t - to)c. 



This shows that the velocity can be determined from the 
coefficient of the term linear in t. 

The Green function has an eigen-function expansion 
for the spatial part 



(16) 



where i?„ are the eigen-values and ipn the corresponding 
eigenfunctions. Since we are considering a dissipative 
system, it is guaranteed that En > 0. For all the eigen- 
functions with En > 0, the Green function will have a 
time contribution G{t,ti, z, zi) ~ exp(— £'„(t — ti)) so 
that in the long time limit, these modes will not con- 
tribute, except for the initial transients. 

The Goldstone mode corresponds to a solution with 
i?o = as can be verified directly with i/jo = d(j)o/dz. 
This zero mode produces the linear term in the solution 
(not an exponential decay in time) and the velocity of 
the front comes from this term only. The velocity then 
is given by 



JMz,){Sf{z,))dz, 
JWoi^iWdz^ 



(17) 
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If the helicase perturbation is operational in a region be- 
, , then the velocity can be written 

AF 



(a) h=0.4d0, W=100 
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(18) 



where £ = J[(f)'Q{z)]'^dz is the kink energy and 

r0+ 



AF: 



<S/(0)d0 



is the free energy cost per unit length in the region of the 
bite. The expression for the kink energy follows from Eqs. 
(fT5)) and ([H]). The velocity found in Eq. [TH]is positive, as 
it should be. 

This procedure can be implemented for the RecG case, 
Eq. [TU) The instability of the interface now leads to a 
moving front but with a negative velocity at least in the 
perturbative regime. The perturbation theory via the 
Goldstone mode yields 



1 

^£ V d 



(19) 



where the limits are the two end points of the bite. What 
we see is a negative velocity because of the reduction of 
the elastic constant and the curvature of the profile of 
the zeroth order interface near (jj = 0. By symmetry one 
would get a positive velocity if the perturbation were at 
6 = 1. 



A. Numerical Solution 

The perturbation theory is around the static solution 
and need not be valid in practical situations. More struc- 
tures can be expected if we look at the nonperturbative 
effects of the helicase term in presence of the dynami- 
cally generated front. To do so, we solved the discretized 
versions of the equations numerically via an implicit pro- 
cedure. The arbitrary lattice spacing Az and time At are 
chosen small for convergence. We have taken At = 10""' 
and Az = 10^^. All the other parameters are chosen in 
these units. 

The results are shown in Fig. |3] and [5] in discretized 
z and time t. The location zq of the interface is taken 
to be the point where cj) = 0.5. We start with a sharp 
interface at time t = 0, with zq chosen away from the 
boundary. Initial transients allow the interface to evolve 
to its natural shape. This shape in the steady state has 
been found not to depend on the initial profile. The 
perturbations are applied over a region as indicated in 
Fig. m For Fig. 2^, the local perturbation favours the 
unzipped state as shown in Fig. [3^ while in Fig. /i^ 
perturbation makes both states unstable (Fig [8)3) . The 
velocity measurements were done once the steady state 
is reached and the interface is not close to the boundaries 
to avoid finite size or boundary effects. 
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FIG. 4: (j) vs. z profile for various times as marked. Both z 
and t axe discretized. (a) For general helicase, Eq. |4]and (b) 
for RecG-type case for Eq. [lOl The horizontal bar indicates 
the location of the perturbation. In (a) the helicase bites at 
(j) = 0.5 while in (b) the bite is at </> = 0.1, in both cases 
over a length 2W. We see the velocity (direction indicated by 
an arrow) towards the zipped side in (a) but in the opposite 
direction in (b). The data are for D = 7 = 1, Hr = 2, so that 
D' < 0. 



A check on the perturbation theory would be the lin- 
earity of the velocity for small perturbations. We have 
verified this and also the fact that for a wide perturbation 
(much larger than the width of the interface) the velocity 
should saturate to the bulk value p^. These verifications 
are not shown here. Fig. 2] shows the change in the direc- 
tion of the velocity by the two perturbations. In (a) the 
helicase bites at a particular value of </> = 0.5, but since 
the perturbation makes the = state mctastable, the 
propagation is always towards the zipped side - opening 
of DNA. In case (b), the bite is also at a particular value 
of (p but close to the boimd state side. We chose (p = 0.1. 
Owing to the instability of the profile, there is a push to- 
wards the unzipped side. The profile in Fig. \^ shows the 
structure developed - the monotonicity of the profile of 
Fig. 2^ is gone. This non-monotonicity can lead to new 
phenomenon not captured by the firstorder perturbation. 
Depending on the width of the bite, there may even be 
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a cancellation in Eq. [19] as we see in the case marked 
400 (ignoring the initial transients) In this situation, the 
helicase would have a stop-and-go type motion, of course 
towards the unzipped side. Though RecG may not do 
so, but there are other helicases like RecQ, RecA which 
are known to have such peculiar motions. Therefore the 
proposed destabilizing mechanism can lead to varieties 
of motions. Moreover, these parameters provide us with 
a set of values to classify or differentiate the helicases 
quantitatively. 




time 



FIG. 5: Velocity as a function of time for the RecG type case, 
because of the continuous change of the profile, the veloc- 
ity may show a time dependence as seen in the case marked 
400 (the spatial width of the perturbation in arbitrary units) . 
Same situation as in Fig. |4]3. 



IV. SUMMARY 

A few comments can be made here. In the case of front 
propagation due to bulk instability, almost a century old 
problem [20|. there is a pushed to pulled type dynamic 
transition as reflected in the nature of the velocity of 



propagation 2l|. In our approach, the biochemically con- 
trollable width W within which a helicase works actually 
acts as a finite length cut-off for the bulk transition, and, 
therefore, it provides a new testing ground for any fi- 
nite size scaling of the bulk nonequilibrium transition. 
This would be an extremely interesting situation where 
a biological problem can shed new lights on an age-old 
physics/material science problem. 

Biological models have so far ignored the importance of 
the interface or treating the Y-fork as a coexistence. We 
hope our results would motivate new experiments to mea- 
sure the elastic energy of the interface (D), possibly via 
independent propagation velocity measurement. In addi- 
tion, the helicases need to be classified by the nature of 
the perturbation and its strength and width. One needs 
to wait for these independent experimental results for any 
quantitative comparisons between the model calculations 
presented here and related biological experiments. 

We summarize our results here. The Y-fork generated 
by a helicase, if kept in equilibrium, is a coexistence of 
two phases. Ordinary helicases work by destabilizing the 
zipped phase near the interface and the local instabil- 
ity leads to the Y-fork motion towards the zipped side, 
thereby opening the fork. For RecG, the aromatic inter- 
action between the helicase and the interface leads to a 
destabilization of the interface, and this could close the 
DNA by a fork motion towards the open side. We show 
this by a perturbation method that also verifies the re- 
sults for the ordinary case. A numerically exact solution 
is used to study the Y-fork motion in presence of the 
dynamically generated front. In the interface destabi- 
lization case, more complex stop-and-go type motion is 
found to be possible. Coming back to the problem of 
replication, our results enable us now to use this formu- 
lation to describe the replication in a co-moving frame 
with the Y-fork, bypassing a direct reference to the heli- 
case. It remains to be seen if one can couple the Y-fork 
motion with the polymerase activity. 
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